Towards a Cold-atom Realization of Quantum Maps with Hofstadter's Butterfly 

Spectrum 
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Quantum systems with Hofstadter's butterfly spectrum are of fundamental interest to many re- 
search areas. Based upon slight modifications of existing cold-atom experiments, a cold-atom re- 
alization of quantum maps with Hofstadter's butterfly spectrum is proposed. Connections and 
differences between our realization and the kicked Harper model are identified. This work also ex- 
poses, for the first time, a simple connection between the kicked Harper model and the kicked rotor 
model, the two paradigms of classical and quantum chaos. 
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The Harper model 1| plays a fundamental role in many 
research areas because it yields the famous Hofstadter's 
butterfly spectrum 0, 0|. This fractal spectrum, first 
discovered in two-dimensional electron systems subject 
to a square lattice potential and a perpendicular mag- 
netic field Q, has found applications in studies of quan- 
tum Hall effect 0, HI, the renormalization group 0, 
high- temperature superconductivity 0], to name a few. 
Various effects on Hofstadter's butterfly spectrum were 
carefully examined [§]. Systems with a butterfly spec- 
trum should be also of general interest to quantum phase 
transition studies because it implies an infinite number 
of phases when some external parameters are scanned. 
Systems with Hofstadter's butterfly spectrum were also 
studied experimentally [1, Q . 

Hofstadter's butterfly can emerge in the quasi-energy 
spectrum of periodically driven systems as well. In this 
context, the kicked Harper model (KHM), adapted from 
the Harper model by considering a delta-kicking poten- 
tial, has attracted vast interests The Hamilto- 
nian of the KHM is given by Ukhm = (L/T) cos(p) + 
K cos(q) S(t — nT), where L and K are two system 
parameters, T is the kicking period, p and q (q S [0, 2w]) 
are conjugate momentum and angle variables, with their 
commutation relation defining the effective Planck con- 
stant h, namely, [q,p] — i% (in a flat phase space, to 
be more precise). The Hilbert space with the periodic 
boundary condition in q is spanned by the eigenfunctions 
\m) of p, with p\m) — mh\m), (q\m) = exp(imq)/\2n, 
and m being an integer. The KHM quantum map as- 
sociated with the unitary evolution for each period T is 
given by 



Ukhm 



-ijt cos(p) -if cos(q) 



(1) 



Because the classical limit of the quantum map Ukhm 
is chaotic in general, this map has become a paradigm 
for understanding (i) how a fractal quasi-energy spec- 



trum affects the quantum dynamics and the associated 
quantum-classical correspondence, and (ii) how the un- 
derlying classical chaos affects the butterfly. For exper- 
imental realizations of the KHM, one early study pro- 
posed to use Fermi-surface electrons in external fields 
[llj . Another study showed that the system of a charged 
particle kicked by a designed field sequence [l2| can be 
mapped onto the KHM. However, these two proposals 
have not led to experiments. Connections between the 
KHM and the so-called kicked harmonic oscillator model 
were also noticed [l3j], but only for the special case of 
K = L. 

Using cold atoms periodically kicked b y a n optical lat- 
tice, about ten laboratories worldwide [3, E, EE E3] 
have realized the so-called kicked-rotor model (KRM) 
[l8| as another quantum map paradigm. Using sim- 
ilar notation as above and in the same Hilbert space 
as the KHM, the Hamiltonian of the KRM is given by 
Hkrm = p 2 /2 + K cos(q) J2 n 5{t — nT). Many variants of 
the KRM, obtained by considering different types of kick- 
ing sequences or additional external potentials, have also 
been achieved. In these studies the experimental setup 
itself has also advanced, from using thermal atoms to us- 
ing a Bose-Einstein condensate (BEC) [lil [lj| that has 
very large coherence width. These ongoing experimental 
efforts motivate the following bold and important ques- 
tion: can these cold-atom laboratories working on the 
KRM also realize the kicked Harper model or its vari- 
ants by slightly modifying their existing apparatus? If 
yes, quantum maps with Hofstadter's butterfly spectrum 
can soon be experimentally realized in many cold-atom 
laboratories, an entirely new generation of experiments 
can be planned, and novel applications of cold-atom re- 
searches may be established. 

Stimulated by our early work seeking a potential con- 
nection between KHM and a variant of KRM [ijj], a 
very positive answer to the above question is indeed pro- 
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vided here. In particular, we show that previous experi- 
mental setup for the so-called double-kicked rotor model 
(DKRM) [I?], H3] already suffices for synthesizing a quan- 
tum map displaying Hofstadter's butterfly, provided that 
one quantum resonance condition therein is met and the 
initial atom cloud is a BEC that has sufficient coher- 
ence width. The butterfly spectrum associated with the 
obtained quantum map is almost indistinguishable from 
the standard result previously calculated for the KHM. 
We then show connections and dramatic dynamical dif- 
ferences between our quantum map and the KHM. In ad- 
dition to experimental interests, the results should also 
motivate more theoretical work on quantum maps with 
Hofstadter's butterfly spectrum. 

Consider then a DKRM that is already experimen- 
tally realized 17] . Using the same notation as above, the 
Hamiltonian of a DKRM can be written as Hdkrm — 
p 2 /2+K± cos(q) J2 n S(t-nT)+K 2 cos(q) 5{t-nT-rf). 
Evidently, in addition to kicks at t — nT, the rotor in a 
DKRM is also subject to kicks at t = nT + n. The asso- 
ciated quantum map Udkrm for a period from nT + 0~ 
to (n + 1)T + Q~ is given by 



- p 2 



-i-T^cos(g)^ (2) 



Udkrm = e^^e"^ cos ^e 

Remarkably, if we now require the parameter T to sat- 
isfy the quantum resonance condition of the KRM, i.e., 
Th = 47r, then due to the discreteness of the momen- 

turn eigenvalues, one obtains e~ %T ^ = 1 when operating 
on any state in the Hilbert space defined above. Under 
this resonance condition we are able to reduce Udkrm 
to Udkrm > 
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FIG. 1: (Color online) Quasi-energy spectrum (denoted e) of 
the quantum map U^ KRM given by Eq. ((3| (top), compared 
with that of the KHM map Ukhm (bottom) given by Eq. fT}. 
Ki/h= K 2 /h = K/h = L/h = l. 



u 



DKRM 



— e lr >* 



cos(g) — M7| 



cos(q) 



(3) 



Because the cold atoms are actually moving in a flat space 
rather than a compact angular space, the quantum reso- 
nance condition is relevant only if the initial quantum 
state is prepared in a definite quasi-momentum state. 
This is certainly within reach of today's experiments. 
For example, two recent experiments |1 6j j studied directed 
transport in a KRM on quantum resonance, with a delo- 
calized BEC (with negligible self-interaction) effectively 
realizing appropriate initial states such as \m = 0). Note 
also that the quantum map U DKRM offers a cold-atom 
realization of a modified kicked-rotor model we recently 
proposed where the kinetic energy term can take 
"negative" values. 

With the U DKRM realized above, we now present one 
key numerical result of this work. Figure [T] displays the 
calculated quasi-energy spectrum of U DKRM by the stan- 
dard diagonalization method [1(3] as a function of fi = t]h, 
for Ki/h — KijTi = I, compared with that of Ukhm as 
a function of %, for K/h = L/K = 1. The map U DKRM is 
seen to generate a beautiful Hofstadter's butterfly. Even 



more dramatically, the butterfly of U L 



DKRM 



resembles the 



previously calculated butterfly of KHM |l(| to such a de- 
gree that the top panel appears to be indistinguishable 
from the bottom panel in Fig. [TJ The generalized fractal 
dimensions, denoted by D 9 , have also been calculated for 
many system parameters, confirming that the spectrum 
is indeed a fractal in general [e.g., for K\ = Ki = 1, 
h = 27r/(l + a), a = (V5 + l)/2, D q=0 « 0.5, same as 
the result for the KHM for h = K, K = L = 1]. Similar 
results have been found for many other system parame- 
ters as well, so long as U = U, K\ = K, and Ki = L. To 
confirm a fractal butterfly spectrum experimentally, one 
may connect the associated characteristics of the quan- 
tum diffusion dynamics (e.g., time dependence of the sur- 
vival probability, the diffusion exponent etc.) with the 
spectrum [21| or attempt to reconstruct the spectrum by 
first reconstructing the time evolving wavefunction. 

Results in Fig. Q] suggest a strong connection be- 
tween a DKRM under quantum resonance and the KHM. 
To uncover this connection let us return to Eq. ([3]) 
and temporarily treat it in a flat phase space without 
the periodic boundary condition. Using the equality 
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FIG. 2: (Color online) Part of quasi-energy spectrum of 
Udkrm in Eq. J3]) and Ukhm in Eq. JTJ for fi = % = 
ho — 267r/41 (the upper two rows) and h = h — 2n — ho (the 
lower two rows). Ki/h = K 2 /h = K/h = L/h = 1. The 
spectral differences between U DK rm and Ukhm are evident. 



f(q + rjp), this treatment leads to 



U 



DKRM 



u 



DKRM 



cos(q+p) g -i-J- cos(q) 



(4) 



where p = r\p is a rescaled momentum variable with 
p = —ihd/dq, Ki — r]Ki, and K 2 — r]K 2 . Equation 
dH) now clearly resembles Ukhm m Eq. |T|), with the 
only difference being that, the first exponential factor of 
Udkrm m Eq- Q contains the cos-function of the an- 
gle plus the momentum, rather than just the momentum. 
This already partially rationalizes the strong resemblance 
between the two panels in Fig. [TJ 

So is the spectrum of U t DK rm identical with that of the 
KHM? Put differently, does there exist a unitary trans- 
formation Q to ensure G'(q+p)G = p, G^qG = <?? If such a 
G exists, then G^U r DKRM G becomes precisely the Ukhm 
in Eq. Q} (with K x -> K , K 2 -> £, and fi^h). Signifi- 
cantly, such a £7 does not exist for the Hilbert space here. 
In particular, the above G transformation is found to as- 
sume the analytical form G — e lq ^ 2h , which violates the 
periodic boundary condition associated with q — > q + 2ir. 
As such, the spectrum of U r DKRM , and hence also the 
spectrum of U 1 DKRM in ([3]), should contain substantial 
elements that are absent in the KHM. 

Indeed, as shown in Fig. 2, a more careful compar- 
ison does expose spectral differences between U r DKRM 
and Ukhm- Motivated by this observation, we are also 
able to find some major differences analytically. For ex- 
ample, for fixed L/h and K/h, the spectrum of Ukhm 
is invariant upon a ?i-change from ho to 2?t — h^. Such 
a symmetry docs not exist in the case of U r DKRM (see 
Fig. 2). A second example is for the special case of 
h = h = Ait. Therein the spectrum range can be easily 
found, which is [-[Kx+Ki)/n, {Ki+k*)/%\ for U r DKRM 
and [—(K + L)/h, {K — L)/K\ for Ukhm (if none of these 
range boundaries exceeds ±7r). 

More insights emerge if we examine one interesting 
classical limit of the DKRM, i.e., the h — > limit (by let- 
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FIG. 3: (Color online) Dynamics of m-variance (denoted a 2 ) 
for quantum maps U 7 DKRM and Ukhm, for the initial state 
|0). The values of h — h are given by ho or 2n — ho, with 
h = 26tt/41. Ki/h = K 2 /h = L/h = K/h = 1. The 
map Ukhm, not U t DK rm, is seen to be invariant upon the 
change ho — * 2-k — ho- Note also that one case associated 
with Udkrm displays localization, and all the others show 
quadratic diffusion. 



ting i] — > with fixed K\ and K%) while keeping Th = An. 
Denote (qi,pi) as a classical trajectory right before the 
Zth kick in the Ti — > limit of U r DKRM . Then one obtains 
P21+1 =J>2i +^ 2 sin(g 2 /); qii+\ = qn +P2(+r; P21+2 = 
P21+1+K1 sin(g 2 i+i); and q 2 i +2 = q2i+\-p2i+2- We stress 
that this classical limit is obtained under quantum reso- 
nance, and is hence unrelated to the direct classical ana- 
log [22j of the quantum DKRM. Upon making a classical 
canonical transformation (q,p + q) — > (Q, P), we obtain 
P21+2 = P21 + K 2 sin(Q 2 ;); Q21+2 = Q21 ~ K\ sin(P 2;+2 ), 
which is precisely the classical map of the kicked Harper 
model. This finding hence firmly binds the butterfly spec- 



trum of U r DKRM with the standard KHM. The emer- 
gence of the classical KHM map from the fi — > limit 
of U 7 DKRM further demonstrates that the spectral dif- 
ferences between U DKRM and Ukhm arise from gen- 
uine quantization effects. Indeed, it is the periodic 
boundary condition in the quantization that disallows 
the above-mentioned unitary transformation G as the 
quantum analogy of the classical canonical transforma- 
tion (q,p + q) -> (Q,P). 

The spectral differences between U r DKRM and Ukhm 
are found to result in profound consequences in the quan- 
tum dynamics. One excellent example is shown in Fig. [31 
demonstrating clearly that the butterfly associated with 
^dkrm violates a symmetry property of the butterfly 
associated with Ukhm- Note also that in one case of 
Udkrm shown in Fig. [3j the quantum diffusion displays 
evident localization that is in clear contrast to the KHM 
dynamics. This localization behavior suggests that the 
width of the sub-bands of the butterfly is effectively zero, 
in agreement with the result shown as the second row in 
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Fig. 2. Indeed, in our numerical analysis the associated 
band width is found to be less than 10 -16 . 

To further motivate interests in the new quantum map 
U r DKRM we also present in Fig. 4 the time-dependence 
of the momentum distribution profile, for a more generic 
h = h that is irrational with it. Results therein show 
again striking differences between U T DKRM and Ukhm, 
especially in that the former case shows a staircase struc- 
ture in the profile. Unlike previous observations of anal- 
ogous staircase structure with a classical origin in an 
off- resonance DKRM [2(J, the staircase structure here 
is purely quantum mechanical. Indeed, both the local- 
ization shown in Fig. [3] and the staircase profile shown 
in Fig. 4 can be related to the unique blocked band 
structure of U 7 DKRM in the momentum representation 
[23l | . The important lesson here is that many impor- 
tant features of a quantum map can be hidden in the 
overall pattern of its butterfly spectrum. For experimen- 
tal interests, we note that one may tune the value of Ti 
and other system parameters to generate different block 
sizes of U r DKRM [23j , thus attaining staircase steps of less 
height and hence more accessible to experiments. For ex- 
ample, for K\ = K2 = 14.4, k ps 1187r/61, the associated 
staircase steps in P{m) have a width of 61, and a height 
2 — 3 orders of magnitude smaller. 

Cold-atom realizations of the non-kicked Harper model 
using static optical lattices were proposed before 
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However, due to the deep lattice approximation therein 
they cannot be extended for the kicked Harper model. 
Based on already available experimental techniques that 
can achieve a double-kicked rotor model tuned on quan- 
tum resonance, here we have proposed a rather simple 
cold-atom realization of a variant of the kicked Harper 
model. The results should open up a new generation 
of cold-atom experiments on quantum maps with a but- 
terfly spectrum. This work also establishes, for the first 
time, a direct connection between the kicked-rotor model 
and the kicked Harper model, arguably the two most im- 
portant paradigms of classical and quantum chaos. 

Acknowledgments: Very constructive comments 
made by Dr. Zai-Qiao Bai are gratefully acknowledged. 
We also thank Prof. C.-H. Lai for his support and en- 
couragement. J.W. acknowledges support from Defence 
Science and Technology Agency (DSTA) of Singapore un- 
der agreement of POD0613356. J.G. is supported by the 
start-up funding (WBS grant No. R-144-050-193-101 and 
No. R-144-050-193-133) and the NUS "YIA" funding 
(WBS grant No. R-144-000-195-123), National Univer- 
sity of Singapore. 



* Electronic address: phygj@nus .edu.sg| 
[1] P.G. Harper, Proc. Phys. Soc. Lond. A 68, 874 (1955). 
[2] D.R. Hofstadter, Phys. Rev. B 14, 2239 (1976). 




10 l 
10-' 
10" 1: 
10" 1 ' 

10- 2 - 
10* 

10~ 31 



: n=4.0xl(f 

'■ n=1.6xl(f jk. 


— n=7.4xl<) 3 
n=2.0xl0' 


- n=3.2xl(f ^ ¥ A'**7 * 




/ // 








1 \ 




J 





-1200 



-600 





m 



600 



1200 
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The staircase structure seen for U^ K rm does not exist for 
Ukhm, despite the fact that their two butterfly spectrum is 
hardly distinguishable. 
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